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Technical Field 



This invention relates to methods and systems 
for maximizing satellite constellation coverage at 
predetermined local times for a set of predetermined 
geographic locations. x 

Background Art 



Space mission design and analysis is typically 
performed to determine constellations that can meet the 
greatest number of mission requirements at the least 
10 possible cost. "During this process, coverageof candi- 
date -constellations is often mat^hjid^ 

Traffic need's^ however, 



nee ds at specif ie d^J^g^t ions - 
could be highly correlated with local time as exempli - 
f i ed^T5y^trhe^aiiy— fee-l-ephon'e"t'rarf"flc wfiTich peaKs^du^ri-ng 
15 t fre— 9L: 0 0 a.m. to 5 : _0_0_p-..m — p e-riod-; Synch-ronization o f 



the^cbnst illation wiETi local time, thus, can increase 
the resource utilization significantly. 



Several existing orbits, such as Sun- synchro- 
nous, orbits, Molniya orbits, etc., do synchronize the 

20 exact passing of the subsatellite point of each individ- 
ual satellite over a predetermined location at certain 
predetermined local times. However, there is no general 
systematic approach to synchronize coverage of an entire 
satellite constellation, which may consist of more than 

25 one satellite, with local time. 
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Disclosure Of The Invention 




It is thus a general object of the present 
invention to provide a method and system for maximizing 
satellite constellation coverage at predetermined local 
5 times for a set of predetermined geographic locations. 

In carrying out the above object and other 
objects, features, and advantages of the present inven- 
tion, a method is provided for maximizing satellite 
constellation coverage at predetermined local times for 

10 a set of predetermined geographic locations. The method 
includes the step of determining a satellite constella- 
tion having a first coverage wherein the constellation 
includes at least one desired satellite having a trajec- 
tory associated therewith. The method also includes the 

15 step of determining a period of rotation for each of the 
desired satellites. Still further, the method includes 
the step of determining a time dependent coverage of the 
satellite constellation based on the period of rotation*, 
and the trajectory of each of the desired satellites 

20 The method further includes the step of tilting the 
trajectory of at least one of the desired satellites to 
obtain a second coverage based on the time dependent 
coverage, wherein the second coverage provides maximum 
coverage at the predetermined local times for the set of 

25 predetermined geographic locations. Finally, the method 
includes the step of generating command signals for 
modifying the trajectory of the at least one desi3re~d~~ 
satellite. 

In further carrying out the above object and 
30 other objects, features, and advantages of the present 
invention, a system is also provided for carrying out 
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the steps of the above described method. The system 
includes a processor operative to determine a period of 
rotation for each of the desired satellites, determine 
a time dependent coverage of the satellite constellation 
5 based on the period of rotation and the trajectory of 
each of the desired satellites, and to tilt the trajec- 
tory of at least one of the desired satellites to obtain 
the second coverage. The system also includes means for 
generating command signals for modifying the trajectory 
10 of the at least one desired satellite. 

The above object and other objects, features 
and advantages of the present invention are readily 
apparent from the following detailed description of the 
best mode for carrying out the invention when taken in 
15 connection with the accompanying drawings. 

Brief Description Of The Drawings 

FIGURE 1 is a diagrammatic representation 
illustrating a satellite communication system of the 
present invention; 

20 FIGURE 2 is a schematic illustration of a 

typical constellation of satellites employed in the 
method and system of the present invention; 

FIGURE 3 is a flow diagram illustrating the 
general sequence of steps associated with the method of 
2 5 the' present invention; 

FIGURE 4 is an equirectangular view of the 
projection of the orbit planes on the surface of the 
earth for a specific satellite constellation; 
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FIGURE 5 is a graph illustrating the number of 
satellites observed per hour at a predetermined geo- 
graphic location before incorporating the method of the 
present invention; 

FIGURE 6 is a graph illustrating the number of 
satellites observed per hour at the predetermined 
geographic location after incorporating the method of 
the present invention; and 

FIGURE 7 is an equirectangular view of the 
projection of the orbit planes on the surface of the 
earth for the specific satellite constellation after 
incorporating the method of the present invention. 

Best Modes For Carrying Out The Invention 

Referring first to Figure 1, a satellite-based 
communications system 10 with a typical geometry for 
practicing the present invention is diagrammatical ly 
illustrated. In general, the system 10- includes a 
plurality of communications satellites 14 belonging to 
a satellite constellation. The satellite constellation 
may consist of any number of orbital planes 15, as shown 
in Figure 2. The satellite constellation of the present 
invention may provide telecommunications coverage, or 
other types of coverage such as, but not limited to, 
remote sensing or surveillance. 

The satellites 14 are in communication with a 
satellite ground station 12, which controls and monitors 
the satellites 14 according to data processed by a 
processor 16. 



• * 



Turning now to Figure 3 , the general steps 
associated with the method of the present invention are 
illustrated. First, a desired satellite constellation 
is determined, as shown at block 20. For example, it 
5 may be desirable to design a satellite constellation 
providing global coverage at all times. Thus, it may be 
desirable to use five satellites in two planes to 
provide such global coverage. It is also desirable to 
have the two planes be more uniformly spaced over the 

10 globe. If the desired geographic location is the border 
of Spain and France, the latitude is 40 'N while the 
longitude equals 0*. The inclination angle of each of 
the satellites may be 45-* . If the desired geographic 
location is Los Angeles, which is near longitude equal 

15 to 120*, the ascending nodes would be 60* and 120° 
correspondingly . 

Each of these satellites are designed to 
follow a specific trajectory defining the orbit of the 
satellite as a function of time. The trajectory of a 
20 satellite orbit is defined by a set of orbital parame- 
ters. These parameters correspond to rotations of the 
trajectory in some inertial coordinate systems. 

Next, the period of rotation of each of the 
satellites is determined, as shown at block 22. If the 
25 point of interest (predetermined geographic location) 
and the orbit are near the equator, i.e., latitude 
equals 0", the orbital period is determined as follows: 

P = [ m s D a D N / (nD N + m s D s ) ] , (Eq. #1) 

where , 

3 0 P is the orbit period with its sign indicating 

whether it is a direct or retrograde orbit; 
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n is an integer with its absolute value equal 
to the number of times that the satellite transverses 
the same geographic longitude within the repeating 
period; 

m s is the number of mean solar day per repeat- 
ing period and must be a positive integer relatively 
prime to n; 

D s is the mean solar day, which is 24 hours or 
144 0 minutes; and 

D N is the nodal day which is the period of the 
earth-rotation relative to the ascending node or any 
point of the orbit plane. If the orbit plane does not 
rotate in the inertial frame, it is the same as the 
sidereal day D x which is about 1436 minutes. 

When the point of interest cannot be covered 
by equatorial orbits, inclined orbits are utilized. The 
orbital period must then satisfy the following con- 
straint : 

T = m 3 D s = m N D N = | (n+m N )P| (Eq. #2) 

where , 

m N is the number of nodal day per repeating 
period which must be a positive integer relatively prime 
to n; and 

T is the repeating period that the coverage 
pattern starts to repeat itself. 

Utilizing the example given above, since the 
desired geographic location is at a latitude equal to 
40* North, Equation 2 must be used to solve for the 
period of rotation of each of the satellites since it is 
not an equatorial orbit. Assuming that the difference 
between m s and m N can be compensated via station keeping, 
n is set equal to 3 and m N = m 3 = 1 . Choosing a repeat- 
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ing period of T equal to 24 hours, P equals 6 hours. 
Thus, the period of rotation of each of the desired 
satellites is equal to 6 hours. 

The geographic longitudinal displacement for 
5 both schemes at the end of the i-th solar day is as 
follows : 

2tt Mod[(i/m a )n,l] . (Eq. #3) 

The geographic longitudinal displacement 
corresponds to the longitude of the subsatellite point 
10 minus the original longitude after the elapsed time 
period. Thus it is guaranteed that the satellite will 
return to the desired longitude at a desired local time. 

Knowing the period of the satellite constella- 
tion, the time dependent coverage provided by the 

15 satellites can then be determined, as shown at block 24, 
utilizing simulation or other similar analysis. Refer- 
ring again to the exemplary satellite constellation, 
which corresponds to the standard ICO (Intermediate 
Circular Orbit) constellation, the constellation con- 

20 sists of two planes of five satellites each at about 
10,355 km altitude and 45* inclination. As discussed 
above, the ascending nodes of the two planes are 180" 
apart at approximately 0" and 180* longitude. The 
f projection of the orbit planes on the surface of the 

25 earth is shown in equirectangular view in Figure 4. 

Based on the period of the ICO constellation, the number 
of satellites observed per hour at the predetermined 
geographic location, e.g., 40° latitude, is shown in 
Figure 5. As can be seen in Figure 5, the satellite 

30 resources are not optimized since a maximum number of 
satellites is seen sporadically throughout the day. 



Since the desired coverage depends on the 
local time at the predetermined geographic location, it 
is desirable to have the maximum possible number of 
satellites providing coverage at the predetermined 
5 geographic location for about 8-12 hours per day during 
the middle of the day. This is achieved by tilting, or 
reorienting, the satellite constellation around the y 
axis in the equatorial plane, as shown at block 26. 
This process is accomplished by repeating the above 

10 simulation after rotating the parameters defining the 
trajectory until the desired coverage is obtained. The 
amount of tilting depends on the constellation. For the 
example given, a tilting of 60° is optimal, as shown in 
Figure 6. The orbital parameters of the tilted constel- 

15 lation corresponds to an inclination angle of 69.3* and 
ascending nodes of 4 0. 9 "E and 13 9.1*E. 

Figures 5 and 6 indicate that the constella- 
tion before tilting has about 5 hours of optimal cover- 
age, namely from 9:00 a.m. to 2:00 p.m. Once the 
20 constellation is tilted, optimal coverage can be 
achieved from 5:00 a.m. to 6:00 p.m., or 13 hours, which 
is 260% better. The projection of the orbit planes in 
equirectangular view after tilting is shown in Figure 7. 

Finally, command signals are generated for 
25 modifying the trajectory based on the desired amount of 
tilting. Since the present invention can be utilized in 
both designing a new satellite constellation or modify- 
ing, a pre-existing satellite constellation, this step is 
accomplished in one of two ways. In the case in which 
30 a new satellite constellation is being designed, the 
tilting operation corresponds to a rotation of the 
original inertial coordinate system. That is, the 



orbital parameters defining the trajectory is first 
expressed in rotation matrices. These rotation matrices 
are then transformed by the desired amount of tilting. 
A new coordinate system is then extracted from the 
tilted rotation matrices to obtain the new orbital 
parameters. The new satellite constellation with the 
optimum coverage is then realized by programming a 
computer 32 (Figure 1) in a launch vehicle 34, such as, 
for example, a rocket or booster, with the new orbital 
parameters. The satellites are then launched into space 
via the launch vehicle with the new orbital parameters 
programmed therein, as shown at block 28. 

For existing satellite constellations, command 
signals must be generated by the satellite ground 
station 12 in order to achieve the desired amount of 
tilting, as shown at block 30. The command signals are 
transmitted by the satellite ground station 12 to 
selected satellites 14 to modify the orbital parameters 
so that the net change of the orbital plane is 60°. The 
command signals instruct the satellites 14 to fire-up 
certain thrusters (not shown) for a certain duration at 
a predetermined time. 

The tilting process of the present invention 
shifts satellite resources away from the equator towards 
the "hot spots" at certain times. Then the motion of 
the earth relative to the orbit plane, not necessarily 
the motion of the earth relative to the individual 
satellite, brings the satellite resources to the daily 
traffic at proper times. Thus, the daily coverage 
provided by the entire satellite constellation matches 
the traffic needs at predetermined local times. The 
overall performance of the satellite constellation is 



# 
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thus improved without any alteration to the space 
segment hardware . 

While the best modes for carrying out the 
invention have been described in detail, those familiar 
with the art to which this invention relates will 
recognize various alternative designs and embodiments 
for practicing the invention as defined by the following 
claims . 



